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REFLECTION OF A WAVE OFF A SURFACE 



BRENDAN GUILFOYLE AND WILHELM KLINGENBERG 



Abstract. Recent advances in twistor theory are applied to geometric optics 
in R^. The general formulae for reflection of a wavefront in a surface are 
derived and in three special cases explicit descriptions are provided: when the 
reflecting surface is a plane, when the incoming wave is a plane and when the 
incoming wave is spherical. In each case particular examples are computed 
exactly and the results plotted to illustrate the outgoing wavefront. 



1. Introduction 



o 

Q 
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C^ ■ In geometric optics, Huygens' principle allows one to describe the propagation 

of light from two alternative perspectives: one can trace the rays or one can trace 
the wavefronts. The drawback with the former description is that following a finite 
number of rays may not describe all of the phenomenon under study, while the 
formation of caustics in wavefronts can cause difficulties in the latter description 

In this paper we utilise recent work [2] in twistor theory to go back and forth 
Cn ' between these two perspectives and thus describe the most elementary of opti- 

cal phenomena: reflection. In particular, we consider the following situation: an 
incoming wave of light is reflected on a surface in M.^ . Given the shape of the incom- 
f~^ . ing wavefront and the reflecting surface, can one describe the reflected wavefront? 

Throughout, we assume that the medium is homogenous and that the speed of 
propagation is unity. 

2 ' The technique we employ in answering this question comes from the minitwistor 
f^ . correspondence which identifies the space of oriented afhne lines in M'^ with the 

tangent bundle to the 2-sphere. This has a long history and has been used in 
various contexts. In particular, it has been used in the construction of minimal 
^\f ' surfaces [8] , solutions to the wave equation [9] and the monopole equation [4] . 

3 . The general context of this work is within the study of line congruences, that 
is, 2-parameter families of oriented lines in M.^ [5] [6] [7]. A line congruence is 
integrable if it is orthogonal to a family of surfaces in R^. At the outset it is not 
clear that the reflection of an integrable congruence in an arbitrary surface is itself 
integrable. However, in Theorem 3 that is precisely what we prove. This is the 
celebrated Theorem of Malus, independently proven by Hamilton in 1827. 

In our formalism, going from an integrable line congruence to the orthogonal 
surface is equivalent to inverting a Cauchy-Riemann operator. The reflection prob- 
lem can then be broken down into four steps. First we transfer from the incoming 
wavefront to its associated ray system and we do the same for the reflecting sur- 
face. The reflection law relates the outgoing rays to these two line congruences at 



Date: July 22, 2002. 



2 BRENDAN GUILFOYLE AND WILHELM KLINGENBERG 

the point of contact, and, finally, we transform the outgoing rays to the associated 
wavefront. 

We present the general formulae relating all of these stages for an arbitrary in- 
coming "wavefront and reflecting surface (Theorem 2). We go on to give explicit 
descriptions for three special cases: when the reflecting surface is a plane (Proposi- 
tion 2), when the incoming wave is a plane (Proposition 3) and when the incoming 
wave is spherical (Proposition 5). In each case, particular examples are computed 
exactly and the results plotted to illustrate the outgoing wavefront. 

The next section contains the background details of the twistor construction 
we use - further details can be found in [3]. In section 3 we deduce the required 
reflection law and establish the fact that the outgoing ray system is integrablc if 
and only if the incoming rays are integrablc. 

In section 4 we turn to the simplest type of reflection, namely, reflection of 
an arbitary wavefront in a plane. Section 5 deals with the general formulae for 
reflection of a plane wavefront in an arbitrary surface. As an illustration of our 
technique we determine the reflection of a plane wave in a sphere and a torus. In 
the final section we give the general formulae for reflection of a spherical wavefront 
in an arbitrary surface, and illustrate the technique in the case where the surface 
is the unit sphere. 

2. The Twistor Correspondence 

The key to our approach is the following geometric construction. Let (x^, a;^, x^) 
be Euclidean coordinates on R^ = C © M and set z = x^ + ix'^, t = x^. 

Consider an oriented line in R^ passing through (z, t) with direction given by 
^ G S*^ C R^, where £, is obtained by stereographic projection from the south pole 
of the unit sphere about the origin onto the plane through the equator. 

Then the minimal distance vector from the line to the origin is given by [2] 









d d 


T^S^, 








r^=}-{z-2ti- 


--zf), 


and the 


point {z, 


. t) is 


a distance 

lz + ^-z+{l 


~mt 



from the point on the line closest to the origin. 

Conversely, given rj and r and the direction of the line ^, the point (z, t) in R^ 
can be found by 

2(7? -^e') + 2^(1 + eC> 



(1+^0^ 



(2.1) 



-2{f^i + riO + {l-et)r 



t= ^'^ '^' 1 ^ ^ ^ . (2.2) 

Now an oriented line in R^ is uniquely determined by ^ and rj as above. This 
is the minitwistor correspondence [4], where the space of oriented lines in R^ is 
identified with the tangent bundle to the unit sphere about the origin. 
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Definition 1. A line congruence is a 2-paramctcr family of oriented lines. 

Definition 2. A line congruence is integrable if it is orthogonal to a family of 

surfaces in M^. 

Suppose the line congruence is given by i^ ^ (CC'^j '^)^ vi^i ^))i ^^^ v G C 
Theorem 1. [3] ^ line congruence is integrable iff 

where d is differentiation with respect to v . 

This equation is the integrability condition for the existence of a real solution r 
to the following equation 

2 {rid^+fjd O 

(i + ea^ 

where r is now considered as a function of v and v . 



g^= T^T;r - (2.4) 



Thus, by the above theorem, given an integrable line congruence, we can find a 
local description for the orthogonal surfaces by inverting the d operator. The real 
constant of integration gives the affine parameter along the lines normal to these 
surfaces. An explicit description of the surfaces in R^ can then be obtained by 
inserting ^, r\ and r, as functions of v and v into (2.1) and (2.2). 

Away from flat points, a surface can be parameterized by its normal direction 
[3], i.e. V = ^ and r\ = F{^,^) for some complex function F. In this case we refer 
to F as the twistor function of the surface, the integrability condition (2.3) reduces 
to the simpler 



_ / 
= a - 


F 


;i + ee)2 


2F 





,(l+eo' 

and the function r satisfies 

dr - 

{i + ^O' 

We will refer to r as the potential function for the surface. 

At a number of points in this paper we use Euclidean motions to simplify the 
equations for reflection. A translation which takes the origin to (z, t) is the qua- 
dratic holomorphic transformation 

e^e, 7j^r] + ^{z-2t^-ze), 

while a rotation about the origin is given by the fractional linear transformation: 



where a, /? G C satisfy aa + [3(3 — 1. 

3. Reflection in a Surface 

Consider a ray with direction ^1 and perpendicular distance vector rji striking a 
surface S at the point (2:0, to)? where the normal direction is ^0 and the perpendic- 
ular distance vector is 770- Let ^2 be the direction of the outgoing wave and 7^2 be 
the perpendicular distance vector (see Figure 1). 
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Normal 



Incoming ray 



Outgoing ray 




Reflection in a surface 

Figure 1 



For the surface S, suppose that the point of reflection is a distance tq from the 
minimal distance point on the normal line to the origin, so that, by (2.1) and (2.2) 



ZO 



2(^0 -%eo) + 2^0(1+^0^0)^0 



to 



-nmio + Vo^o) + {l-^Uo)ro 



(1 + CoU^ (l+eo^o)^ 

On the other hand, since both the incoming and outgoing rays contain (zo, to), 

iOT i= 1,2. 

Combining the previous three equations we have, after some rearrangement, 

(1 + U.?VO - (eo - ^^?fjo + (Co - 6)(1 + ?"o6)(l + eoCo)ro 



Vz = 



We turn now to the law of reflection. 



(l+Co^o)^ 



(3.1) 



Proposition 1. If S,i is the direction of the incoming ray, 1^2 the outgoing ray 
direction and ^0 the normal direction at the point of reflection, then 
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2eoCi + 1 - Co^o 
(l-Coeo)Ci-2eo' 



(3.2) 



Proof. A rotation about any point ^0 on P^ is described by a unitary fractional 
linear transformation: 



Let S^g be the rotation 



Rot5o(0 

? == 1. Th 

Rot^-;(0 

5?o(0 = 



P^ + a 
where a, /3 G C satisfy ad + /3/3 = 1. The inverse rotation is given by 



f3^ + a 



(3.3) 
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where 



Clearly S^^ is a rotation that takes ^o to zero (the North pole). Thus, if we denote 
reflection through ^o by R^g, then R^^ = S7^ o Rq o S^^. 

Now -Ro(C) — ~Cj so using (3.3) and the definition of S^^ we find that 

„ ... ^ (Cogo - l)e + 2^0 

Alternatively, we can define i?^„ as above and check that it is a rotation (obvious), 
and satisfies R^o{(,a) — Co and R^g o R^g{(,) = S,- Thus it is a rotation through 180" 
about ^Q i.e. reflection in ^o- 

Finally, the law of reflection says that the direction ^2 of the outgoing ray is 
obtained by reflecting the antipodal direction of the incoming ray S,i through the 
normal direction ^ at the point of reflection. 

To complete the proposition, we note that the antipodal map on P^ is ^ -^ — C^^; 

so that 

, „ f l\ 2^0^1 + 1-^0^0 
« - R^o 



V (l-6eo)ei-2eo' 
as claimed. D 

The equations governing reflection are: 

Theorem 2. Consider a wavefront given by vi — > (Ci(z^i, Pi), ?7i(z^i, j>i)) reflecting 
off a surface given by i>q -^ (Co(^o, i^o), ^o(^'Oi '^o)) 'J'^c? ro(t'o, f'o)- Then the reflected 
wavefront is determined by 

/- _ 2CoCi + 1 - CoCo /o A\ 

^^"(1-^0^0)6-2^0' ^'-'^ 



m 



m 



(1 + Coeo'r/o - (go - iiYm + (Co - a)(i + eo6)(i + eoeo)r-o, 

(i + eoeo)2 

(Co - iifm - (1 + i^lifm + (Co - a)(i + Coei)(i + eoeo)ro 

((1-^0^0)6-2^0)2 



(3.5) 



(3.6) 

Equation (3.5) determines the intersection of the incident rays with the surface, 
while (3.4) and (3.6) determine the direction and perpendicular distance from the 
origin of the reflected rays. 

Proof. These come from combining the reflection law (3.2) with (3.1), after some 
rearrangement for 772. D 

In general, the reflected line congruence can be parameterized by vi, the param- 
eter for the incoming wavefront. 

Theorem 3. The surface, incident and reflected congruences satisfy the following 
relationship: 



•^1 ~T-, — , A ;r N/-, — , ^ r N '^O 



(1 + 66)2 (1 + 66)2 'V (i + 66)(i + 66) 



(3.7) 
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where di is differentiation with respect to vi, the parameter for the incoming wave. 
Proof. We start by differentiating the reflection law (3.4): 

-(1 + eoCo)'aiCi + 2(e"o - h?diio + 2(1 + ioh?diio 



ai6 
ai6 



[(1 - ^0^0)6 - 2Co]2 

-(1 + eoeo)'aiei + 2(^"o - ei)'ai6 + 2(1 + ioiifdiio 



[(1 - eo^o)a - 2^0]^ 

The reflection law (3.4) also implies 

(i+eoeo)'(i+a6) 



1 + 66 = 



[(1 - 66)6 -2eo][(l- 66)6 -26] 
Thus, using (3.5) and (3.6) we find 

V2di^2 + fi2di^2 ^ vidi^i + fjiBi^i 2(|6-6P- |i + 66n(^o^i6 + ??o^i6) 
(1 + 66)2 " (1 + 66)2 (1 + 66)3(1 + 66) 

, 2[(6 - 6)(i + 66)[(i + 66)5i6 - (i + 66)ai6]] 



(1 + 66)2(1 + 66) 
2[(6 - 6)(i + 66)[(i + 66)5i6 - (i + 66)ai6]] 



ro 



^0- 



(1 + 66)2(1 + 66) (3.8) 

Now, since we are reflecting in a surface, we have the potential function tq satisfying 

2(?7o9o6 + ^o9o6) g 

(1 + 66) 

and so 

2(770^16 + ^0^16) _ 2(?7o(5ii^o5o6 + ^it^ot^o6) + rjo{diVodo£,o + BiVoBo^o) 
(1 + 66)2 - (1 + ^^^-^)2 

_ 2(9iPo(?7o9o6 + ^oc^o6) + diiyo^Tjado^o + fjodg^g)) 

(1 + 66)2 

= diPgdgro + 9ii/o9oro 
= dirg. 
Substituting this in (3.8) we get that 

mdi^2 + ^2^16 ^ ??i'9i6 +??i^i6 I6-6|2 - |i + 66|2 k 
(1 + 66)2 " (1 + 66)2 (i + 66)(i + 66) '''° 

, 2[(6 - 6)(i + 66)[(i + 66)^16 - (1 + 66)5i6]^ 



(1 + 66)2(1 + 66) 
2[(6 - 6)(i + 66)[(i + 66)^16 - (1 + 66)9i6] 



''o 



^0, 



(1 + 66)2(1 + 66) 

and the last three terms on the right hand side are di of the real function, as stated 
in the theorem. D 

As a corollary we get the Theorem of Malus: 

Corollary 1. A reflected congruence is integrable if and only if the initial congru- 
ence is integrable. 
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Proof. Since the second term on the right-hand side of (3.7) is di of a real function, 
the integrabihty condition (2.3) means that the outgoing congruence is integrable 
if and only if the incoming wave is integrable. D 

In the next sections we consider three special cases: that of reflection of an ar- 
bitrary wavefront in a plane, and that of plane and spherical wavefronts reflected 
off an arbitrary surface. In all cases the resulting wavefront can be explicitly de- 
termined. 

More generally, assume that the incident wave can be described by rji = -Fi(^i, ^i) 
for some complex function Fi. Suppose further that we can invert (3.5) for ^i as a 
function of ^o ^nd 770- Then we can substitute this in (3.4) and (3.6) to find both 
^2 and ?72 as a function of ^0 and tjq. Finally, we can find the potential function r2 
by integrating (2.4). 

4. Reflection in a Plane 

Consider an arbitrary wave described parametrically by ^i(i^, 9) and rii{i>, v). We 
want to determine the resulting wave after reflection in a plane. By a rotation we 
can align the plane with the x^x^-plane so that the initial wave lies in the region 
x^ > 0. 

Proposition 2. The reflection of a wave given by ^i(j^, iy) and r]i{h', v) in the x^x^- 
plane is 

6 = i, m = -|- (4.1) 

Proof. The x^x^-plane is given by ^0 = ''0 = 0- Inserting this in the reflection law 
(3.4) immediately gives the first of the above equations. Inserting it in (3.5) and 
(3.6) we get 

m = m- S.fm, m^m- -^- (4-2) 

SI 

Adding the first of these to ^1 times its complex conjugate gives 

_rti+ &rji 

Inserting this in the second of (4.2) gives the result. D 

4.1. Reflection of a Spherical Wavefront in a Plane. Consider a spherical 
wavefront with source (0,0, ii). This congruence is given by 771 = — 2ii^i. Substi- 
tuting this in (4.1) we find 

, 1 2*1 

Combining these two we see that 772 — 2ii^2, which is a spherical congruence with 
source (0,0, — ti). Thus we retrieve the well-known law of reflection, whereby the 
reflection of a spherical wave centered at (0, 0, ii) in a plane mirror is another 
spherical wave with virtual centre (0, 0, — ti). 
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5. Reflection of a Plane Wavefront 

Consider an incoming wavefront with fixed direction, i.e. a plane wavefront. In 
this case, the reflected wave can be expressed explicitly in terms of the reflecting 
surface: 

Proposition 3. The reflection of a plane wavefront with direction ^i off a surface 
given by ^o(^: ^) o-'^^d '?o(^7 ^) 's given by 

^'"(1-^0^0)6-2^0' ^'-'^ 

_ (go - $i)'^o - (1 + Ui?m + (go - Ci)(l + 6ei)(l + io^o)r, 
"^^ ((1-66)6-26)2 ' (5 2) 

the reflected congruence being parameterized by the parameter value v of the surface 
at the point of reflection. 

Proof Since the incident angle is constant, there is no need to determine the points 
of intersection, and the above equations are just the equation of reflection (3.4) and 
equation (3.6). D 

In this proposition, the resulting wavefront is parameterized by the point of 
reflection on the surface. For a plane wave it is also possible to parameterise the 
outgoing wave by its direction: 

Proposition 4. Suppose S a surface in M.'^ , given by the rjQ = -Fo(6 C) with potential 
function 7*0(667 ^■^ struck by a plane wave with normal direction — ^j~ . Then the 
reflected wave is given by the line congruence (S^, 772 = -^2(66/' with 

^2(66 = 7 ((1 + 7)'J^o(6,6) - (e - 76)'^o(6,6) + 2(6 - 67 ro(6,6)) , 

4 (5.3) 

where 



, _ 666 - 1 ± v/(i + 6)(i + 66)(i + g6)(i +16) .... 

^° 6i + 6i) + 6(i + C6) ' ^ ' 



and 



^ = J (l + ^^"i)(l+^^") 



(i + C6)(i + 66) 

The two signs yield the same line congruence with opposite orientation, and can be 
chosen to point in the outgoing direction. 

Proof. From the reflection law (3.4) with 6 = C ^-nd 6 ^ — Cr ^^ have 

, ^ (66 - 1)6 + 26 
266 + 1 - 66 ■ 

This can be recast as a quadratic equation for 6j which has solution (5.4). From 
these, after some calculation, we find the following expressions 



1+66 - ± ^^ , 1-66 ^^ , 

1 + 60 = ^ , 6-e = ^ , 
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where we have introduced 

6 = 1- e^'CiCi /? - ^(1 + ai) + 6(1 + eCi)- 

Then we compute the following combinations of the above 



i + gg"o _ ±V(i + e6)(i + eO + v/(i + e6)(i + aei) 

1 + ^0^0 2v/(i + e6)(i + aei) 



Co-? _±6V(i + eei)(i+e?)-?V(i + ea)(i+aei) 



i + CoCo 2v/(i + ea)(i + a6) 



(eo-e)(i+eco) ^^i(^^_^^^/ (i+ea)(i+^e) 



i + eoCo 2^^^ V (i + ^?i)(i+^i^i) 

Substituting these in (3.6) gives the result stated. D 

We also have the following corollary for waves travelling down the x'^-axis: 

Corollary 2. Suppose S a surface in MP, given by the rjQ ~ -Fo(CiC) with potential 
function ro(a^), is struck by a plane wave moving down the x'^-axis. Then the 
reflected wave is given by the line congruence (^, ri2 — F2{S,,£,)) with 

F2{t O^UU + V^l + Ctj ^o(Co, Co) - em^Ao) - 2^\fl+i^ ro(Co, Co) J , 

where 



, -1 + Vl + gg 
Co- J . 

Proof. This follows from the above Proposition 4 by setting ^ = 0. D 

By way of example, we now compute the twistor function for a plane wave after 
reflection off the unit sphere and the rotationally symmetric torus. 

5.1. Reflection of a Plane Wavefront in the Unit Sphere. By translation 
we can move the unit sphere to the origin in M'^ and by a rotation we can set the 
direction of the plane wave to be down the x^-axis. Thus we can use corollary 2. 

The twistor function and potential for the unit sphere at the origin are _fo = 
and ro = 1. Thus, substituting these in (5.5), the reflected wavefront is given by 



F2 = -UJl + H, r2 = J-^ + C. 

Figure 2 shows the resulting outgoing wavefronts (which are axially symmetric, but 
not spherically symmetric). 
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5.2. Reflection of a Plane Wave front in a Torus. Consider the torus with core 
radius a and meridian radius b which is rotationally symmetric about the x^-axis. 
This has twistor function and potential 



^-|(^"Wf 



To = = + 0. 



In what follows we will consider the case where a — 2 and b — 1. We will deal with 
three cases: a wave travelling down the x'^-axis, up the x^-axis and one striking the 
torus at an angle of 45° to the vertical. 



5.2.1. Plane wavefront moving down the x^-axis. The line congruence resulting 
from reflection of a plane wave moving down the x'^-axis can be found directly 
from corollary 2. The result is 



Fo 



2(1 - CO - v^v/TTI? R 



r2 






where we have integrated F2 to obtain r2. 

Figure 3 shows the complete wavefront as it leaves the torus. Notice the annular 
shadow. 




Plane wave reflected in the torus 

Figure 3 
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5.2.2. Plane wavefront moving up the x^-axis. In order to determine the reflection 
of a plane wave approaching from another direction we shaU use Proposition 3. If 
the direction of the incoming wave is ^i , substitution of the torus twistor function 
and potential in (5.2) gives 



m 



m - a)'eo - (i + CoCO'CoKi - CoCo) + (Co - 6)(i + eoei)[(i + Coeo) va^+ 4CoC"o] 



VC^((l-eo^"o)Ci- 2^0)2 
Integrating this gives the potential function 



(5.6) 



7'2 



4[2(ieo - Cip - |i + eoeinv/eo^+ +(eoeo(i - aa) - ^6 - eoei)(i + eoCo)] 



(i + ei6)(i + eoeo)2 



(5.7) 



For a wave travelling along the x^-axis we set ^i = 1. The resulting wavefront 
is shown in Figure 4. 




Plane wave along the x-axis 
reflected off a torus 



5.2.3. Plane wavefront at an angle of 4^ to the vertical. Finally, a portion of the 
reflected wavefront generated by an incident wave making an angle of 45'^ with the 
cc^-axis is shown in Figure 5. This is obtained by setting ^i — 2.4 in (5.6) and (5.7). 




1 2 



3 4 



Plane wave at an angle of 45° 
reflected off a torus 
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6. Reflection of a Spherical Wavefront 

Assume that the wavefront is spherical, that is, the rays have a single focus. By 
a translation we can move this focus to the origin and then the wavefront is simply 
given by r/i = 0. Every point in M.^, other than the origin, is contained on two 
oriented lines in the congruence (the same line with both orientations). 

Proposition 5. A spherical wavefront, with focus the origin, reflected off a surface 
(not containing the origin) determined by S,oii^,i^), riQ{v,i>) and ro{iy,P) gives rise 
to the congruence 

, ^ 2r?o + 2^0^^ ± 2go/3o 

^' 2(^0^0 - Cor/o) - (1 + Co^o^ro ± (1 - eoeo)/3o ' ^ ' ' 



whe 



ai 



a2 



m ^ aim - a2'?o - 2(1 + 6Co)q:3?'o, (6.2) 

2f?oCo - (1 + ioio)ro ± /3o 



2(6?7o 


- Co%) 


-(l+eoCo)Vo±(l- 


-CoCo)/3o' 




2r/o + 


?o(l+eoCo)'^o±eo/3o 




2(^o??o 


- iom) 


-(l+eoCo)2ro±(l- 


-CoCo)/3o' 


4?7o?7o?o - 


- im - J?oCo)(l + ^o^oVo ± {vo 


+ %eo)/3o 



"' (2(^0% - ?"or/o) - (1 + ^0^0)^0 ± (1 - ?oeo)/3o)2 ' 
and 

Po^\/4.mm + {l + io^o)^rl 

Here the two different signs give the same line congruence with the opposite orien- 
tation, and can be chosen to point in the outgoing direction. 

Proof. Setting 771 = in (3.5) we get a quadratic equation for ^1, with solution: 
2(%eo + Uo) - (1 - ^Ui)ro ± (1 + eoeo)/3o 
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2(^0 -%eg) + 2^0(1 + ^0^0)^-0 



Inserting this in (3.4) yields the expression for ^2 in the theorem. In addition we 
find that 

(27^0^0 - (1 + ^o^o)ro ± /3o)(l + Co^o) 



a -6 







1 + Coa = 



2(?7o-J7oeo) + 2^0(1 + eoCo)'-o ' 
(2ryo + eo(l + ^o^o)ro ± eo/3o)(l + ^^o) 



while 



(1 - ^0^0)6 - 2^0 = 



2(%-??oC§) + 2^0(1 +aC"o)ro 
(2(Co^o - eo'?o) - (1 + ^0^0)^0 ± (1 - CoCo)/3o)(l + Coeo) 



2(??o- ^0^0) + 2^0(1 + ^0^0)^0 
Substituting these in (3.6) completes the theorem. D 
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6.1. Reflection of a Spherical Wavefront in the Unit Sphere. Consider the 
unit sphere, centred at (0,0,-2). This is determined by 



'70 = 2^0, 



ro = 1 - 2 



1 - ^0^0 
1 + CoCo ■ 



Thus 



and by (6.1) 



6 



2^0(2(1 + eoeo)+/3o) 



l-2^oeo-3egCg + (l-eoeo)/3o' 
where we have chosen the + sign in the equations. Further straightforward, if 
lengthy, calculations establish that 



ai 



1 + Co^o + /3o 



a2 = 



_ go(3(l + eoCo)+/3o) 
l-2CoCo-3eo'^"2+/3o' "' l-2^ofo-3^2^"2+/3o' 

U^+Cdo + M 
"' (l-2eo?o-3eo2^"2+/3o)2- 
Finally, substitution of these in (6.2) gives 



m 



4^0(1 - 3eoeo)(l + 3eo?"o + Vl + lOeoCo + ^Qeo) 



1 + ^0^0 " 7^o'Co + 9CoCi? + (1 - mo + s^ieo)VT+im^^+W^ 



To find the potential function r2 we invert (2.4) with v = S^q, since we are parame- 
terising the outgoing rays by the direction of the normal to the surface at the point 
of intersection. After some computation we obtain 



r2 



-2(1 - 3^oeo)Vl + 10^oeo+9?2Co' _^ ^ 



1 + 11^0^0 + 19^0^^1+9^^1 

At this juncture we have found parametric expressions for ^2, V2 and r2, that 
is, the complete description of the outgoing wave. If we graph the wavefronts the 
results are similar to the plane wave reflected off the sphere. By way of comparison. 
Figure 6 compares a cross-section of the two wavefronts in a plane containing the 
a;^-axis. The shaded area represents the shadow cast by the surface. 

Reflections in a Unit Spliere 





Plane Wave 



Spherical Wave 
Source (0,0) 
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